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Abstract 

We point out that the non-gaussianity arising from cubic self interactions of the inflaton field is 
proportional to £iV e where £ ~ V" and N e is the number of e-foldings from horizon exit till the 
end of inflation. For scales of interest N e = 60, and for models of inflation such as new inflation, 
natural inflation and running mass inflation £ is large compared to the slow roll parameter e ~ V' 2 . 
Therefore the contribution from self interactions should not be outrightly ignored while retaining 
other terms in the non-gaussianity parameter /nl- But the N e dependent term seems to imply 
the growth of non-gaussianities outside the horizon. Therefore we briefly discuss the issue of the 
constancy of correlations of the curvature perturbation £ outside the horizon. We then calculate 
the 3-point function of the inflaton fluctuations using the canonical formalism and further obtain 
the 3-point function of We find that the N e dependent contribution to /nl from self interactions 
of the inflaton field is cancelled by contributions from other terms associated with non-linearities 
in cosmological perturbation theory. 
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Non-gaussianities arising from self interactions of the inflaton field are usually treated 
as negligible compared to those arising from non-linearities in cosmological perturbation 
theory [1, 2] because they are proportional to higher order slow roll parameters, such as 
£. In this paper we point out earlier results indicating that this contribution includes a 
term proportional to N e , the number of e-foldings after the scale of interest has left the 
horizon, which for our horizon scale is approximately 60 by the end of inflation. Then for 
certain models of inflation with relatively larger higher order slow roll parameters, such as 
new inflation, small field natural inflation (/ < 1.5Mpi), and running mass inflation, this 
contribution can be comparable to other contributions from non-gaussianities in cosmological 
perturbation theory. 

However the dependence on N e would seem to imply that n-point (n > 2) correlations of 
the curvature perturbation ( grow outside the horizon. Hence we discuss the related issue 
of constancy of ( outside the horizon. The derivations regarding the constancy of (k are 
classical. They only imply that the 2-point function of ( to lowest order, ~ |0c| 2 <5 3 (k — k'), 
is constant outside the horizon. (Hatted quantities denote quantum operators.) For higher 
point functions of (, there is an additional requirement associated with convergence of an 
integral over time, as we discuss below. 

Having clarified this we then work in a gauge in which 5(f), the fluctuation in the inflaton 
field, is not zero and recalculate the 3-point function of 5<f)k, using the canonical formalism 
which we believe has not been done before. This agrees with results of earlier calculations 
using the path integral approach and field equations. 

We then relate 54>k to (k- We argue that the quantisation of the relation between ( and 
5<p derived from the SN formalism is not ideal if we wish to study any possible growth of 
n-point functions of £ outside the horizon. Using a slightly different relation which explicitly 
includes dependence on the final time, we then obtain the 3-point function of (k, including 
the term associated with self interactions of the inflaton field which is proportional to N e . 
Our calculation of (C 3 ) is straightforward and does not involve complicated field redefinitions. 
Finally we find that the time dependent self interaction contribution is cancelled by other 
terms associated with non-gaussianities in cosmological perturbation theory. 
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For an inflaton with a cubic interaction the bispectrum parameter /nl is 
6 



5 /NL 



^ + 7 - N e + ^3 |^ ^ fcifcj - 

where N e is the number of e-foldings of inflation from the time the mode of interest leaves the 
horizon at t cx till the time t, which can be at any later time during inflation. Our expression 
is similar to that in Eq. (38) of Ref. [3] with their t* replaced by t, and we have replaced 
N* by —N e = —H(t — t cx ). We shall explain the distinction between our expressions later. 
The slow variation in H can be ignored in N e since f^L is to first order in slow roll, e, r\ 
and £ are the slow roll parameters evaluated at t. 

1 [V'V 1 2 



(2) 

2 V V J 2 H 2 K ' 

r] = — ~ — ^ + e 3 
' V H(j) 

We use £ rather than £ 2 as used by some authors and M P = M P \/\/8ir has been set to I. 
7 0.577216 is Euler's constant. 

The contribution due to the cubic interaction of the inflaton is the £ dependent term 
above and was obtained in Refs. [3-6]. The argument to ignore this term is that it is 
proportional to £ ~ V" and is hence negligible compared to terms proportional to e ~ V' 2 
and i] ~ V". This is valid for models such as chaotic inflation for which e > rj > £. 
However, in general, the hierarchy is rj > £ > o > while e may be larger or smaller 
than other slow roll parameters. In particular, for new inflation, small field natural inflation 
(/ < 1.5Mpi), and running mass inflation, £ > e. In general, for small field models with a 
concave-downward potential, e <, 0.0001 [7]. Therefore the £ term should not be outrightly 
ignored while retaining other terms in /nl above. 

There is another reason why the £ term in /nl is not automatically small compared to 
other terms. For scales of the order of our horizon today N e is about 60 if t is at the end of 
inflation. Therefore the term £N e in / NL above can be large. For a new inflation potential 
of the form V = V - fi(f) 3 ((f) > 0), £ = 0.5 rf. Also, r] = 0.5(n s - 1) + 3e w 0.5(n s - 1) [8]. If 
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we take n s = 0.96 [9] then r\ is -0.02 and £iV e is 0.012. e is much smaller. Thus we see that 
£iV e is comparable to rj and e <C £N e . 1 2 

The presence of the time dependent N e in / NL seems to contradict the notion that n- 
point functions of ( do not grow outside the horizon. It may be argued that since the 
curvature perturbation ( does not grow outside the horizon the N e contribution from the 
non-gaussianities of the inflaton can not contribute to /nl- Note, however, that this term 
has been obtained independently in Refs. [3-6]. So we now consider the literature on the 
constancy of ( in the context of non-linear cosmological perturbation theory. 

Salopek and Bond [11] first introduced the generalisation of the Bardeen-Steinhardt- 
Turner variable ( [12] for the case when one wishes to consider non-gaussianities in the 
curvature perturbation. The Salopek-Bond ((x) is constant outside the horizon. In other 
works, such as Refs. [13-18], the constancy of ( is shown for a (k mode. The above works 
deal with classical ( and the constancy of classical ((x) or (k only implies that the quantum 
2-point function (CkiCk 2 ) = (27r) 3 |Cfc 1 | 2 <^(ki — k 2 ) is constant outside the horizon (at lowest 
order). 3 But it is not obvious that other higher n-point functions of ( are constant outside 
the horizon. 

One may argue that after horizon exit the curvature perturbations are classical and so 
the constancy of ( should imply constancy of n-point function for n > 2 also. As argued 
by Weinberg in Ref. [19], the perturbations are classical in the sense that commutators 
involving ( and its time derivative go to for large t. However, this implies that quantum 
n-point functions of zeta can nevertheless grow outside the horizon, but only as powers of 
In a and not as powers of a. It is the time dependence due to In a ~ iV e for the 3-point 
function which is the focus of this article. 

1 Interestingly, Appendix B of Ref. [10] obtains a similar result but erroneously concludes that the N e 
dependent result of Ref. [4] and the N e independent result of Ref. [1] are the same because they are of 
the same order. 

2 In Ref. [3] it is argued that evaluating expectation values at the end of inflation may not be valid for 
large N e ?a 60 because of divergences of the form e m+2 7V™ (m > 1). However for the potential we are 
considering e is much smaller than 1/60. 

3 Higher order corrections from loops can give time dependent or N e dependent contributions as mentioned 
in Sec. VI of Ref. [19]. 
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The n-point functions of ( are given by [1, 19] 

00 ft rtN 

(6(t)) = J2 iN dtN dtN-i--- / dh 

N=0 *° *° *° 



X 



[^ J (* 1 ),[^ J (t 2 ),..-[^ / (t JV ),6 / (t)] ■■■]]) , (5) 



where O(t) can be any product of £ operators, Oi(t) is 0{t) in the interaction picture 
generated by the quadratic part of the Hamiltonian, and t is some early time. Note that 
the expectation values are obtained in the in-in formalism and so the bra-s and kets refer to 
j n (0| and \0) in respectively. Hj is the interaction Hamiltonian and includes terms that are 
third or higher order in (. For the three-point function at lowest order this reduces to 

(C\t))= i J*dt>([H I (t%Cht)]) (6) 

If (k is constant, ( ~ e lk ' x CfcCfc + e^'^Q^k ls constant. But for (C 3 (^)) to be constant outside 
the horizon one must ensure that the contribution to the integral above from to t is 
suppressed. Note that ( is related to 5<f) and ((50) ) grows outside the horizon. 

The convergence of the integral for large t, and certain other conditions for the constancy 
of (( n ) outside the horizon have been discussed in general in Ref . [20] . But in Ref . [20] (see 
Eq. (29)) only gaussian fluctuations of the inflaton are considered. 4 So one should verify 
whether or not (( n ) is indeed constant outside the horizon when one includes non-gaussian 
fluctuations of the inflaton. We now check this explicitly for the three point function of 
( while including a cubic interaction of the inflaton field. The 3-point function of ( has 
been obtained by other authors. Largely, the self-interactions of the inflaton are ignored. 
Moreover, assuming that there is no contribution outside the horzion, the integral in Eq. 
(6) is cut off at t*. In cases where one first relates £ to 5<p using the 5N formalism and 
then calculates the 3-point function, obtaining any contribution from evolution outside the 
horizon is precluded by the adopted formalism. 

3 

We work in a gauge in which 5<p 7^ 0. We first calculate ((50) ) using the equivalent of 
Eq. (6) for the 3-point function of 8(j). Our results agrees with those obtained in Ref. [3] 
using field equations. We then relate £ to 5(f) and use this to obtain (<^ 3 ) and /nl- The N e 
dependent term mentioned above associated with the cubic self interaction of the inflaton 



4 Note that Ref. [11] also considers only gaussian fluctuations of the inflaton, since 5(j> is set equal to H/ (2tt) 
in the evaluation of ( in Sees. IIIB and HID. 
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does appear in /nl- The 3-point function (( 3 ) has also been obtained in a gauge in which 
8(f) is [1, 21]. Since these calculations do not include self-interactions of the inflaton their 
results do not include the N e term. 

The relevant part of the action S for 5<p can be expressed as the sum of terms quadratic 
and cubic in 5<p. For notational convenience we hereafter replace 5<p with Q, and let 
represent the background homogeneous field. Then 



S — So + S3 



(7) 



where [22, 23] 

So = / dt d 3 x a 3 



Q 2 



(8) 



and £3 is given in Eq. (3.6) of Ref. [1]. Retaining terms to leading order in slow roll 
parameters we obtain 



S3 



dt d 3 x a 3 



AH 



QQ 



\2 



\4jQ(diQ) 2 - ~ 2 d t Qd^Q - V'(0)Q 3 
a 2 AH a 2 6 



where 



d 2 ^ 



2H 



0Q 



(9) 



(10) 



In Ref. [1] the last term in S3 above, which is proportional to £, was ignored. As we have 
argued earlier the contribution of this term could actually be larger than that of other terms 
above for certain models of inflation. Therefore it ought not to be ignored at this juncture. 

The shift function JVj = d{ip, and ip differs from \ of Ref. [1] by a factor of a 2 . Using eq. 
(10), S3 can be rewritten as 5 



S3 



dt d 3 x a 6 



AH 



One can obtain if, and thus Hj, from the lagrangian in S 2 + S 3 . We follow Ref. [24] 
for dealing with the Q dependent interaction terms. After obtaining H T (Q, Uq) we replace 
Q and Uq by Q in and U in respectively, and then set ITj n = a 3 Qi n . Keeping terms upto 



5 The action £3 provided in Ref. [22] may contain typographical errors. The prefactor for the last term of 
S3 in Eq. (53) of Ref. [22] should be a~ 2 rather than a -4 , and the definition of ip and in Eq. (43) is 
not in agreement with eq. (2.24) of Ref. [1]. However, Eq. (54) for ip is correct. 
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first order in <p/H, as in the action, we then get Hj(Q,Q) = —L int . Hereafter we drop the 
subscript in. 

The contribution to (Q(ki,t)Q(k 2 ,t)Q(k 3 ,t)) from each term in Hi is given below. 
We use the conformal time r, defined by dt = adr, instead of t and let the initial time 
correspond to r = — oo. The final time corresponds to the reheat time. 



1. The QQ 2 term 



h = H)Qki(r)Q k 2(r)Qk3(r)(2n) 3 5 3 (J2^ 

i 

+ C.C. 



dr' 



dr 1 



(12) 



Q k is given by [25] 



Qk = i 1 ~ exp ('5? 



'1 + kr) exp (—ikr) 



iH 



(13) 



kV2k 

which reduces to iH/(ky/2k) for \kr\ <C 1. 6 There are 6 permutations of k\,k2,ks for the 
expression within the integral. Then 



dr'd)(r' 



k\k\ (1 - ik X T'y ktT ' + perm + c.c. (14) 



where kt = k\ + ki + k%. Replacing the lower limit, — oo, by — oo(l — i5) and setting 5 to 
after taking the limit eliminates the contribution of the lower limit. Integrating by parts, 
taking the limit k(r <C 1, and using the complex conjugate to avoid listing some terms, we 
get 

% HHt) 



h = - 



4n 



/Co/bo . 

' " [ ik 



r) ik^jj) ihrfjj) <t>"{r) 



(ik f Y 



Uk f ) 2 



Sihfir) ihrf^r) 0""(r) 

/ • 7 \ A / • 7 \ A 



(iktY 



{iktY 



dr'\-htk l( t> m \r') + (1 - ik 1 T , )cf>""\T , )]}- : - \ + perm 



+ c.c. 



(15) 



6 There are higher order corrections to the late time mode functions, as discussed in Ref. [21]. 
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where 0' = d<p/dr and so on. 

In the Appendix we assess the higher derivative terms. We find that the term is 
proportional to rfe 2Ne . Similarly the term is proportional to rfe me . These higher order 
terms in slow roll parameters are (increasingly) larger than the terms proportional to (for 
t] = 0.02 and N e = 60). However it is not consistent to consider them here as we have 
ignored terms higher order in slow roll parameters in the action. But this is an indication 
that there may be convergence issues at higher orders in the slow roll parameters and these 
will have to be handled with care. 7 Similar behaviour may be expected while working in 
the 5(f) = gauge as, for example, in Ref. [21] where integrals for (£ 3 ) include powers of 
e ~ O.5 2 /if 2 in the integrand. 

Having noted our concern above, we hereafter do not include terms higher order in slow 
roll parameters. Then 



/l = - 2X II^ (27r)353 E^ (r) 



2 3 " / l" / 2 3 



2JU2 



kt 



+ 



+ perm 



(16) 



The prefactor of 2 comes from the complex conjugate. There are a total of 6 permutations 
of the variables (k±, k 2 , k 3 ). The interchange of k 2 and k 3 gives the same expression as above. 

2. The Q{diQf term 
h = (-i)Qki(r)Q k 2(r)Q k3 (r)(2n) 3 5 3 (J2^ 



f 

j — ( 



dr'a 2 ^[(-k 2 ■ h) QUt'M^QUt') + perm] + c.c. 
i H 3 (t) 



-k 2 ■ h 



dr'(j)(T' 



)(1 - ijbir')(l - ik 2 r')(l - ik 3 T')e iktT ' + perm 



+ c.c. 



^ kikj kik 2 k 3 

~ k t + / ^ ~i — i — n — i- P erm 

k t kf 



(17) 



Clearly there is a symmetry in (k 2: k 3 ) interchange. k 2 ■ k 3 can be replaced by (k\ — k\ — kf)/2 

— * 

using h = 0. Note that the first term above is obtained from the real part of % exp[ik t r]/r 



7 A concern regarding using perturbation theory in slow roll parameters may also be found in Ref. [3], as 
mentioned earlier. 
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in the limit kfT <C 1. 



3. The diQid^Q)^ term 

h = ^)Qki(r)Q k2 (r)Qk3(r)(2n) 3 6 3 (J2^) 



f 



dr i a 2 <t> h-k 2[ ^ ( ^dQ* k2 (T r )dQi z (T r ) 



2H k\ 



= 2- 



1 H s < 

( 

1 H 3 (t) 



= 4n^) (arW ?^ T) 



dr' 
ki ■ k 2 

k\ 



h perm] + c.c. 



dr 

pp kik 2 k 2 
2 3 " / l" / 2 3 



+ 



+ perm 



(27r) 3 5 3 (^^)0(r)(r i -fc 2 ) 



/to /Ci /Co 

!T + -^+ perm 



(18) 



The 3-point function of Q is ii + I 2 + h plus the contribution from the cubic self- 
interaction. For the self-interaction contribution we use the expression given in Ref. [3] 
which agrees with Refs. [4, 6]. Using Mathematica to include all the permutations and then 
simplify their sum gives 



(Q(k 1: t)Q(k 2: t)Q(k 3: t)) 



(2tt) 3 5(A; 1 + k 2 + k 3 ) 



H 2 V" 



4IM 



( ~l k t +k t J2k i k j + l{ 1 - + 1 + \n | k t r\ } k lj 

\ i<j i / 



+ 



# 4 



1 / 1 



= (2TrY5(k 1 + k 2 + k 3 ) 



H 2 V" 



^^^ k i k j Z tfkjkk 



(19) 



x 



+ 8TPF^ ( 2 ^ " h ^ ^ " 2 ^ 



(20) 



where i, j = 1,2, 3, k t = ^ ki, and we take all ki have approximately the same magnitude. 
t = —l/(aH). The term proportional to V", which is due to the cubic self-interaction of 
the inflaton, was obtained in Refs. [3-6]. In \k t r\ = ln[(aH) cx / (aH)] f=s m[a ex /a] = — N e , and 
thus one gets an N e dependent term. 8 This term is not explicitly cancelled by any other 



The result in Ref. [5] differs by a sign and a term [3]. In addition, In | Ait 77 1 in Ref. [5] should be set to 
— N e and not +N e . Ref. [26] also obtains a V" term with a In a dependence. 
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term in the expression above. Moreover, as we show later, this term is also not cancelled by 
the time variation of slow roll parameters in other terms. 

The above uses the action/Lagrangian and the canonical formalism to calculate the 3- 
point function of Q. Ref. [22] uses the path integral formalism to obtain the 3-point function. 
(It does not consider the contribution from the cubic self interaction term.) Ref. [3] uses 
the solutions of the Heisenberg field equations to obtain (Q 3 ), including the self interaction 
contribution, and our results agree with Eqs. (20) and (29) of Ref. [3]. Ref. [27] explicitly 
shows the equivalence of the form of (Q 3 ) obtained in Refs. [3, 22]. 

In the SN formalism, the gauge invariant quantity t) is the difference in the number 
of e-foldings of evolution between some time t* and t at x and the number of e-foldings 
between t* and t for an isotropic homogeneous background, where t* lies on a spatially flat 
slice of spacetime with field values (f>(x,t*) while t belongs to a spacetime slice of uniform 
energy density, is typically chosen to be a few e-foldings after the relevant scale has left 
the horizon. 

C(x,t) = N[p(t),<j>(a!,Q]-N\p(t),<i>(U)] 
dN[p(t),<f>(x,Q] 



+ 



d(f)(x, U) 
1 d 2 N[p(t),(f)(x,Q} 



5(f)(x, t* 

4>{u) 



2 d(f)(x,uy 



5<j)(x,Q 2 + ■■■ (21) 

4>{U) 



where (f>(t*) is the spatial average value of at t*, 6<f>(x,t*) = <f>(x,t*) — (f>(t*) = Q(x,t*), 
and • • • refers to higher order terms that have been omitted. (We have temporarily rein- 
troduced (f)(x,t)). Dependence of iV on is ignored in the slow roll approximation 
N[p(t),(f)(x,U)] is given by 

N\p(t),</>(x,U)] = [ H[(f)(x,t)}dt 



r mm w$ (22) 



Then 



<9iV[p(t),0(f, Q] 
d(f)(x, t*) 



H[<t>{U)] 
<t>{u) <f>(t*) 
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d 2 N[p{t)^{x,U)} 



d(j)(x, t* 



d (V{<f>) 



«u) ~ d( t> \ V '^) 



-1-^ (24) 



The + (— ) in the first equation is for V'((f>) > (< 0), or equivalently, for < (> 0). 
Below we will consider the sign as for V'(4>) < as relevant for new inflation. (However the 
various contributions to (( 3 ) are ultimately dependent on only even powers of 0, and so the 
final expression for (£ 3 ), or / NL , is independent of the sign of V'(4>).) Then ( is related to 
Qby 

C(k, t) = ±^Q(k, U) + \ (l - |r) J Q(h - q, U)Q(q, **) + •••, (25) 

The + (-) in front of the first term is for V'{4>) > (< 0), or equivalently for <p < (> 0). 

In the 5N formalism, ( is independent of t, and there is no t dependence on the r.h.s. 
of Eq. (25). If one directly quantises the above relation, as in Refs. [3, 22], then Eq. (25) 
implies 

(CA,t)cA,t)CA,t)) 

= -- — -3{Q(h,QQ(k 2 ,QQ(k 3 ,Q) 
(2e*)2 

/ ^Q&-q,QQ(q,Q+perm) (26) 

The r.h.s. above does not depend on t but that is because of the way ((t) was defined as in 
terms of Q(t*). 

One may instead argue that the quantum field should be expressed as a function of 
quantum operators at t. Furthermore, t and t* are defined on different hypersurfaces and 
so correspond to different 'definitions' of time. We need a relation between the quantum 
operators ((t) and Q(t) which is valid at all times t, including prior to and after horizon 
exit, and both operators should be functions of the same t. For this one may use Eq. (A8) 
of Ref. [1] and 



C(k, t) = ±^=Q& t) + \ (l - |) j ^ Q(h - q, t)Q(q, t) + -.., (27) 

'• • • ' includes terms quadratic in Q but which will not contribute to the three point function 
because their coefficients are suppressed for t > t ex , i.e., outside the horizon. (As discussed 
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in Sec. 3 of Ref. [1], ( and 5(f) are defined in different gauges, i.e, they are functions of 
different time variables t (uniform density gauge) and t (spatially flat or uniform curvature 
gauge) respectively. But t — t + T(t, x) and so Q on the r.h.s. of the equation above can be 
expressed in terms of the same time variable as on the l.h.s.) Then 

(C(4*,)C&,*)C(*3,0> 

= -—^(Q(ki,t)Q(k 2 ,t)Q(k 3 ,t)) 
(2e)3 

+~(l-^)(Q(ki,t)Q(k2,t) J ^Q(h-q,t)Q(q,t)+ perm) (28) 

The first term can be evaluated using Eq. (20). The expectation value in the second term 
above is 

h = (2tt) 3 5 3 (^ h) [2Q kl (t)Q* kl (t)Q k2 (t)Ql 2 (t) + perm] (29) 

i 

= 8e\27r) 3 5 3 (J2 kMCWCfo) + perm] (30) 

i 

The factor of 2 on the first line above comes from different ways of contracting the Qs, and 
there are 3 permutations involving ki,k 2 ,ks. P^ is defined as 

m - (3D 

The terms dropped in Eq. (27) would have contributed similar to the second term in Eq. 
(28), i.e., without any time integral as for the first term, and hence would be evaluated only 
at a time t > t ex , when their contribution will be suppressed as mentioned above. Defining 
/nl as in Ref. [3] 

(C(*i,f)C(*2,f)C(*3,f)> = (27r) 3 5(fc 1 + fc 2 + fc 3 )^NL^P C (A;,)Pc(^), (32) 

i<j 

one gets / NL as in Eq. (f). (/ NL as defined in eq. (32) is the negative of / NL in Refs. [1, 22].) 
The key difference between our expression for / NL and that in Ref. [3] and other works is 
that our /nl is a function of t rather than t* because we expressed ((t) as a function of Q(t) 
instead of Q(t*)- 

Our calculation of (( 3 ) above has been rather straightforward. We neither use extensive 
integration by parts to rewrite the action nor do we invoke any field redefinition as in Refs. 
[1, 21, 22]. Furthermore, if t corresponds to the end of inflation then N e ps 60 for modes 
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entering the horizon today. Then, as mentioned earlier, the £N e term in /nl is comparable 
to the rj term and a priori should not be ignored. 

Now we consider the constancy of the 3-point function of (. To determine how / NL 
changes during inflation after a mode crosses the horizon, we take the derivative of /nl with 
respect to t. Now using d/dt = <pd/d(p we get 

^ ~ [4e 2 - 2 V e] H 

^ ~ [4e£ -tf-<T]H (33) 

where a = V' 2 V""/V 3 . Then df Nh /dt « (5/6) d[-$,N e -r)]/dt = 0. Here we have kept terms 
to first order in slow roll parameters as in /nl- Thus for the terms considered above df^/dt 
is zero indicating that /nl(^) = /nl(4x), i-e., the N e factor in £iV e , which corresponds to 
growth outside the horizon, is cancelled by changes in other terms and does not contribute 
to the bispectrum. 9 However note that (Q 3 ) in Eq. (20) is a function of terms proportional 
to £ and e and so the N e dependent term in (Q 3 ) is not cancelled by the time variation of 
other terms. 

The argument above for the constancy of £ closely follows that in Ref. [3]. In Ref. [3] it is 
argued that /nl is constant outside the horizon by taking the derivative of /nl(£*) effectively 
with respect to t*. One may argue that to study evolution of (C 3 (^)) outside the horizon one 
should take derivatives with respect to t. But then df^/dt will be by construct, as in the 
SN formalism ((k,t) is expressed in terms of 5<f>(k,t*) and /nl is independent of t. 

In conclusion, we have calculated the 3-point function of the curvature perturbation ( 
using the in-in formalism in a gauge in which 5<p ^ 0. The calculation of the 3-point function 
of ( has generally been done ignoring inflaton self interactions, and using the 5N formalism 
which by construct does not allow one to check for evolution outside the horizon. We have 
included the contribution associated with non-gaussian fluctuations of the inflaton due to 
a cubic self interaction, and this is proportional to £iV e and grows outside the horizon. If 
we take N e to be the number of e-foldings of inflation after the mode of interest has left 
the horizon till the end of inflation, then for our current horizon scale N e is 60. In new 



9 Note that these arguments are for adiabatic fluctuations of the inflaton and do not apply to any non- 
adiabatic fluctuations or fluctuations of other scalar fields considered in Refs. [4-6]. 
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inflation, small field natural inflation and running mass models of inflation, e < £ < rj, and 
£iV e is then comparable to other contributions to the non-gaussianity parameter /nl, and 
should not be outrightly ignored. The N e dependent term corresponds to evolution outside 
the horizon. However on including the time dependence of other contributions to /nl this 
time dependent growth cancels. Our results also indicate that there may be issues related 
to the convergence at higher orders in perturbation theory. Higher order contributions can 
render /nl to have stronger time dependence. 
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Appendix 

In this Appendix we compare the 0(r), 0(r)',0(r)", 0(r) w and 0(r) w/ terms in Eq. (15). 
We will ignore terms proportional to e and derivatives of e which are small for inflation 
models that are of our interest. We will also use Eqs. (2,3,33). Then rj ps — 0/(iJ0), 
H = —eH 2 ps 0, and rj ps —£H. We take fcj ~ k t . 

The terms within curly brackets in Eq. (15) are ~ <j>/kt- 

The 0' term is ~ 



k\T ■ , ,, k\T ki 6 1 - 

-TT^) = "TT^ = -TSTF « J<h ■ ( 34 ) 
ft, t h, t h t n h t 

Thus this is smaller than the terms. 

The <t>" terms are ~ <p"/kf = a 2 '<p/kf. Now ps -Hr](j). Therefore ps £<j)H 2 + rf(j>H 2 . 

Then the 0" terms are ~ 

^H 2 + rj 2 H 2 )j> (35) 



„2 



, 2 , ,2 ; / !\ 



^H 2 + rj 2 H 2 



a 2 H 2 v y K 

= exp(2iV e )(£ + r ? 2 )^ (36) 

where we use fc t = a ex H cx and ignore the variation in H. With £ = 0.5 rf, rj = —0.02 and 
N e = 60, the above is approximately 10 49 times larger than the term. 
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The <p'" term is 



kiT u„ hr g--;- ha 2 ---- , 

1£* =J4 a + = -Hi4* (37) 

Now w (£ + r] 2 )H 2 <j). Therefore ps — (4£ + r] 2 )r]H 3 <j). Then the 0"' term is ~ 

-^-(4^ + r? 2 )?^ = e 2 ^(4£ + r? 2 )ri (38) 



which is a factor of 77 less than the 0" term. 

The 0"" terms are ~ W/ /A; t 5 . Using 0' from above, ^ ^ (ll^ 2 + 4£ 2 + r] 4 )H 4 <j). Then 
the 0"" terms are ~ 

^-10^^ = 10e 4 N^ (39) 
fcf dt ai x H 4 x k t 1 k t y ' 

For rj = —0.02 and N e = 60, the above term is a factor of 10 98 larger than the term. It is 

also larger than the 0" term by a factor of 10 if exp (2N e ). 

Thus terms higher in slow roll parameters in evaluating the integral in Ii are larger than 

the lowest order terms. This also holds for I 2 and J 3 . 
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